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Report  Abstract 


The  subject  program,  conducted  through  tight  coupling  between  theory,  experiment  and 
computation  and  reported  in  20  journal  articles,  has  focused  on  the  chemistry  and  dynamics  of 
laminar  and  turbulent  flames  of  surrogate  jet  fuels  in  environments  simulating  various 
operational  aspects  of  aero-engines.  For  studies  on  combustion  chemistry,  we  have 
experimentally  acquired  such  chemistry-affected  data  as  the  ignition  criteria  and  laminar  flame 
speeds  of  fuel-air  mixtures,  which  are  respectively  relevant  for  low-  and  high-temperature 
chemistries,  and  used  these  data  to  assist  the  development  of  their  reaction  mechanisms. 
Auxiliary  contributions  include  the  formulation  of  a  fitting  formula  for  the  falloff  curves  of 
unimolecular  reactions,  which  is  more  accurate  and  easy  to  evaluate  than  the  classical  Troe 
formula,  and  the  development  of  a  Chemical  Explosive  Mode  Analysis  (CEMA)  computation 
algorithm  that  allows  on-the-fly  assessment  of  the  local  reactivity  of  a  mixture  and  thereby 
identification  of  the  dominant  modes  of  burning  and  flame  stabilization  dynamics. 

For  studies  on  flame  dynamics,  we  have  identified  and  classified  various  modes  of  intrinsic 
flame  front  instability  for  laminar  flames,  including  thermal-diffusive  cellular  and  pulsating 
instabilities,  as  well  as  the  (cellular)  hydrodynamic  instability  which  is  promoted  at  elevated 
pressures.  The  potential  of  self-turbulization  through  flame  front  wrinkling  has  also  been 
assessed.  For  turbulent  flames  we  have  developed  a  novel  apparatus  allowing  the  study  of 
expanding  spherical  flames  in  controlled  turbulence  and  at  constant,  elevated  pressures.  The 
measured  turbulent  flame  speeds  show  accelerative  propagation  and  correlation  with  an 
appropriately-defined  turbulent  Reynolds  number  based  on  global  flame  radius  and  Markstein 
diffusivity. 
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Research  Summary 


PART  I:  DEVELOPMENT  OF  COMPREHENSIVE  CHEMICAL  KINETIC  MECHANISMS  FOR 

LARGE-SCALE  COMPUTATIONS 


1.  Premixed  and  Non-premixed  Ignition 

The  interest  in  the  study  of  nonpremixed  ignition  of  hydrocarbons  is  to  understand  the 
controlling  mechanisms  and  responses  of  ignition  in  a  strained  flow  field,  in  which  chemistry 
and  transport  phenomena  are  inherently  coupled.  The  parameters  governing  ignition  not  only  are 
of  practical  interest  but,  together  with  those  in  other  systems,  such  as  shock  tubes,  flow  reactors, 
well-stirred  reactors  and  laminar  flames,  form  an  extensive  matrix  for  the  development  and 
validation  of  comprehensive  hydrocarbon  oxidation  mechanisms.  In  the  past  few  years,  we  have 
been  systematically  and  extensively  acquiring  experimental  data  and  conducting  computational 
simulation  of  the  ignition  of  a  number  of  fuels  in  the  laminar  counterflow,  with  pressure  up  to  a 
few  atmospheres.  While  the  extent  of  fuel  variation  investigated  is  quite  extensive,  it  is 
nevertheless  still  rather  insufficient  in  that  they  are  primarily  gaseous  under  atmospheric 
conditions.  To  extend  the  investigation  to  fuels  that  are  liquid  under  room  conditions,  a  new 
experimental  set  up  for  the  nonpremixed  ignition  in  stagnation  flow  was  designed  for 
experimentation  at  elevated  and  reduced  pressures,  ranging  from  0.61  to  3.0  atm..,  and  with 
quantified  flow  fields,  respectively  recognizing  the  importance  of  pressure  on  the  ignition 
kinetics  and  the  need  to  use  the  local  strain  rate  to  quantify  the  characteristic  flow  time. 
Measurement  of  the  axial  velocity  profile  was  conducted  with  LDV  over  the  local  strain  rates  of 
60  to  350/s.  As  in  the  counterflow  ignition  of  gaseous  hydrocarbons,  the  ignition  temperature  of 
n-heptane  increases  monotonically  with  increasing  strain  rate  and  decreasing  system  pressure. 

Numerical  simulation  of  the  ignition  response  with  detailed  and  reduced  chemistry  and 
transport  was  also  conducted.  The  mole  fraction  of  n-heptane  at  the  liquid  surface  was  found  to 
be  insensitive  to  the  variations  of  strain  rate  and  pressure  examined  in  the  present  study,  and  its 
slight  variation  was  considered  to  have  minimal  effect  on  the  ignition  response.  A  detailed  high- 
temperature  oxidation  model,Mech-130,  yielded  better  agreement  with  experiments  than  a 
smaller  model,  Mech-58,  at  0.61  and  1  atmosphere,  while  the  comparison  is  unclear  at  1.5  and 
3.0  atmospheres  [Fig.  1-1].  An  explosive  mode  analysis  was  conducted  using  the  kinetic 
mechanism  of  Mech-130  to  facilitate  understanding  of  the  ignition  mechanism.  A  transition  from 
radical  explosion  to  thermal  explosion  was  clearly  demonstrated  by  the  explosion  pointers  [Fig. 
1-2].  An  explosion  participation  index  identified  the  controlling  reactions  leading  to  ignition. 
Reactions  of  the  H2/CO  chemistry  and  those  involved  C2-C3  radicals  were  found  to  be  important 
at  the  ignition  states. 
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Figure  II- 1:  (a)  Measurement  and  simulation  of  ignition  temperature  versus  strain  rate;  (a) 
p=0.61  atm ,  (b)  1.0  atm ,  (c)  1.5  atm  and  (d)  3  atm. 
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Figure  1-2:  Explosion  pointer  versus  oxidizer  temperature;  p=latm,  Vboundary  =50  cm/s, 
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Following  the  work  on  n-heptane,  we  conducted  a  study  to  acquire  experimental  data  on  ignition 
temperatures  in  nonpremixed  diffusive  system  and  laminar  flame  speeds  of  CH4/C2H4/N2 
mixtures.  These  data,  together  with  numerical  simulations  using  selected  kinetic  schemes, 
facilitate  our  understanding  of  the  low-  to  intermediate-temperature  ignition  chemistry  and  high- 
temperature  flame  chemistry  of  ethylene-enriched  methane,  respectively.  Computational  analysis 
was  then  performed  to  identify  the  controlling  elementary  reactions  in  ignition  and  examine  how 
the  oxidation  pathways  shift  with  the  relative  fuel  concentrations.  Results  show  that  ethylene 
addition  enhances  the  ignitibility  and  flame  propagation  of  methane.  Furthermore,  both  the 
ignition  temperature,  Tign,  and  the  laminar  flame  speed  respond  steadily  to  variations  of  the 
ethylene  concentration  [Fig.  1-3].  The  experimental  measurements  were  then  compared  with  the 
simulation  results  to  scrutinize  the  low-  and  high-temperature  chemistry  of  the  mechanisms, 
associated  with  ignition  and  flame  speeds  respectively.  Results  show  that  both  mechanisms  yield 
largely  similar  results  for  both  the  ignition  temperature  and  laminar  flame  speed  [Fig.  1-3]. 


5 


Mixing  Ratio,  a 

Figure  1-3:  Ignition  temperature  of  CH4/C2H4  mixture  versus  mole  fractions  of  ethylene,  (a) 
p=latm,  k=300/s.  (b)  p=5atm,  k=300/s ,  (c)  Laminar  flame  speeds  of  CH4/C2H4  mixture 
versus  mixing  ratios  at  various  system  pressures,  cp=1.0 


The  effects  of  thermal  stratification  (through  both  mean  and  fluctuations  in  initial  temperature) 
and  turbulent  mixing  timescales  on  the  auto-ignition  of  a  lean  homogeneous  n-heptane/air 
mixture  at  constant  volume  and  high  pressure  were  investigated  by  direct  numerical  simulations 
with  a  new  58-species  reduced  n-heptane/air  kinetic  mechanism.  In  the  first  parametric  study,  the 
homogeneous  ignition  delay  was  held  constant,  and  twelve  cases  with  varying  initial  mean 
temperature  straddling  the  NTC  regime  were  studied  with  different  degrees  of  temperature 
fluctuations  imposed.  The  displacement  speed,  Damkohler  number,  and  chemical  explosive 
mode  analyses  verify  that,  in  general,  larger  T’  induces  greater  temporal  spreading  of  the  mean 
heat  release  rate  because  the  deflagration  mode  is  predominant  at  the  reaction  fronts  for  large  T’. 
However,  spontaneous  ignition  prevails  for  small  T,  and  hence,  simultaneous  auto-ignition 
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occurs  throughout  the  entire  domain,  resulting  in  an  excessive  rate  of  pressure  rise.  For  mean 
temperatures  lower  than  the  NTC  regime,  e.g.  with  To  =  850  K,  the  ignition  delay  is  increased 
with  increasing  T’.  On  the  contrary,  the  ignition  delay  is  significantly  reduced  with  T’  for  mean 
temperatures  greater  than  the  NTC  region  (e.g.  To  =  1008  K).  For  mean  temperatures  within  the 
NTC  region,  e.g.  with  T0  =  934  K,  the  combined  effects  of  high  and  low  temperatures  manifest 
themselves  such  that  the  ignition  delay  is  increased  for  small  T’  but  is  advanced  with  large  T’ 
[Fig.  1-4].  In  the  second  parametric  study,  the  homogeneous  ignition  delays  were  halved  and 
doubled  for  cases  with  To  =  1067  K  and  To  =  754  K,  respectively.  For  the  cases  with  high  To  = 
1067  K,  the  ignition  delay  is  significantly  reduced  and  the  temporal  spreading  of  the  mean  heat 
release  rate  is  also  enhanced  with  increasing  T,  similar  to  the  cases  with  To  =  1008  K.  For  cases 
with  low  To  =  754  K,  however,  the  ignition  delay  and  temporal  spreading  of  the  mean  heat 
release  rate  are  not  significantly  affected  by  thermal  stratification.  This  is  because  the  turbulence 
timescale  (=2.49  ms)  is  short  relative  to  the  second-stage  ignition  delay  time,  and  hence,  thermal 
fluctuations  are  homogenized  by  fast  mixing  prior  to  the  second-stage  ignition.  Therefore,  the 
second-stage  ignition  is  not  significantly  affected  by  turbulence. 

Finally,  in  the  third  parametric  study  the  ratio  of  the  turbulence  to  ignition  delay  timescale  was 
also  found  to  change  the  ignition  characteristics  of  the  mixtures.  A  fast  turbulence  timescale  is 
able  to  homogenize  the  mixture  such  that  ignition  is  more  apt  to  occur  by  spontaneous  ignition 
for  both  large  and  small  T  cases.  Longer  turbulence  timescales,  however,  are  not  able  to 
homogenize  temperature  fluctuations  prior  to  thermal  runaway.  Therefore,  for  cases  with  small 
T\  ignition  occurs  independently  in  each  ignition  kernel,  and  for  cases  with  large  T\  ignition 
occurs  primarily  by  deflagration,  resulting  in  a  smooth  mean  heat  release  rate.  Overall,  the  effect 
of  different  turbulence  timescales  on  the  ignition  delay  is  small  compared  with  that  of  thermal 
stratification. 
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Figure  11-4:  Homogeneous  ignition  delay  at  constant  volume  with  initial  pressure  of  40  atm 
as  a  function  of  initial  temperature. 

2.  Studies  on  Chemical  Kinetics 

Unimolecular  and  the  reverse  recombination  reactions  play  an  important  role  in  many  chemical 
reaction  systems  related  to  atmospheric  chemistry  and  combustion.  Evaluation  of  these  reaction 
rates  is  however  computationally  demanding  because  of  their  complicated  pressure  and 
temperature  dependence,  especially  for  large  reaction  systems.  Approximate  formulas  of  high 
accuracy  are  therefore  desired  in  the  simulation  of  these  systems.  Our  recent  work  [Zhang  and 
Law  (2009)]  approximated  the  fall-off  curves  of  unimolecular  reaction  rate  constants  by 
proposing  a  fitting  formula  based  on  the  Kassel  integral,  which  is  the  unimolecular  reaction  rate 
expression  of  the  classical  RRK  theory.  Compared  with  previous  formulas  such  as  that  of  Troe 
[Troe  (1979),  Gilber  et  al.  (1983)]  and  Prezhdo  (1995)],  this  (ZL)  formula  not  only  is  more 
accurate  in  the  correlation,  but  it  also  has  a  simpler  mathematical  form  that  is  computationally 
more  efficient  to  evaluate.  In  this  period  of  the  project,  we  have  extended  the  formulation  by 
accounting  for  the  effects  of  tunneling  on  unimolecular  reaction  rates  and  the  corresponding 
falloff  curves.  A  fitting  formula  for  the  pressure-  and  temperature-dependence  of  unimolecular 
reactions  with  tunneling  is  subsequently  proposed.  Compared  to  the  Troe,  Prezhdo  and  ZL 
formulas,  the  present  formula  shows  substantially  improved  fitting  precision  in  the  tested 
reactions  C2H3  — ►  C2H2+H,  C2H5  — >  C2H4+H  and  LC4H3  — >  C4H2+H,  for  which  tunneling 
becomes  significant  at  different  pressures  and  temperatures.  As  an  improved  version  of  the  ZL 
formula,  the  present  expression  includes  only  one  additional  fitting  parameter  to  describe 
tunneling.  Furthermore,  the  tunneling  correction  factor  also  has  a  mathematically  simple  form  for 
computational  expediency. 
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During  this  period  we  have  also  focused  on  hydrogen  abstraction  reaction  mechanisms.  The 
hydrogen  abstraction  reaction  of  alkenes  by  the  hydroxyl  radical  is  dominant  at  temperatures 
above  700  K,  because  of  the  concomitant  rapid  reduction  in  the  rate  of  the  addition  reaction  due 
to  fewer  and  less  efficient  stabilizing  collisions  of  the  adducts  which  tend  to  decompose 
thermally  and  rapidly  before  further  reaction.  The  H-abstraction  reactions  of  butene  isomers  by 
the  OH  radical  were  investigated  by  quantum  mechanical  and  kinetics  approaches.  A 
comparative  study  of  the  geometry  optimization  and  vibrational  frequencies  of  the  transition 
state  structures  of  s-allylic  H-abstraction  was  performed  and  it  was  found  that  results  of  the 
BHandHLYP/6-311G(d,p)  method  is  comparable  to  those  from  the  CCSD(T)/6-31G(d)  method 
with  efficient  computational  time.  Consequently,  the  BHandHLYP/6-311G(d,p)  method  was 
applied  to  determine  the  electronic  structures  of  the  stationary  points  of  the  potential  energy 
surface  and  the  minimum-energy  paths  (MEP)  of  H-abstraction  reactions  of  isomeric  butenes  by 
the  OH  radical.  To  increase  the  accuracy  of  the  computed  potential  energy  surfaces,  single-point 
energy  calculations  were  applied  by  using  the  CCSD(T)/6-31 1++G(d,p)  method. 

3.  Flame  Stabilization  and  Chemical  Explosive  Mode  Analysis  (CEMA) 

Flame  stabilization  is  essential  in  the  understanding  of  lifted  jet  flames.  Depending  on  the  flame 
configuration,  various  stabilization  mechanisms  have  been  proposed  based  on  premixed  and  non- 
premixed  flames,  auto-ignition,  and  turbulence-flame  interactions.  Some  postulations  involved 
are,  however,  difficult  to  verify  with  experiments,  particularly  when  there  is  limited  access  to 
detailed  spatial  and  temporal  species  information.  Numerical  simulation,  particularly  direct 
numerical  simulation  (DNS),  provides  abundant  information  associated  with  the  complex  flow 
and  species  fields.  We  have  used  DNS  data  for  lifted  turbulent  jet  flames  set  to  extract  the  salient 
information  on  the  controlling  stabilization  mechanism  from  the  massive  datasets. 

We  have  developed  the  chemical  explosive  modes  in  chemically  reacting  flows,  based  on  the 
concept  of  CSP,  as  a  new  diagnostic  tool  to  delineate  ignition  and  flame  fronts  in  complex  flows. 
An  abrupt  change  in  the  mixture  explosivity  was  employed  for  the  detection  of  premixed  flame 
fronts,  which  are  difficult  to  detect  using  conventional  methods  involving  only  individual 
parameters  such  as  temperature,  mixture  fraction,  heat  release  rate  and  species  concentrations. 

A  singular  eigenvector  matrix,  or  defective  chemical  Jacobian,  was  observed  at  the  crossover 
point  of  the  chemical  explosive  mode.  It  was  found  that  this  singularity  was  caused  by  the 
convergence  of  the  explosive  mode  to  the  conservation  modes.  A  species  explosion  index  was 
consequently  defined  to  normalize  the  diverging  explosion  pointers  near  the  singularity.  Radical 
explosion  and  thermal  runaway  were  then  distinguished  with  the  explosion  pointers  for  auto¬ 
ignition.  The  present  method  was  first  tested  on  homogeneous  auto-ignition  and  one-dimensional 
laminar  premixed  flames  and  subsequently  applied  to  analyze  the  DNS  results  of  a  turbulent 
lifted  hydrogen  jet  flame  in  a  heated  coflow.  Lean  laminar  and  rich  turbulent  flame  fronts  were 
detected,  and  the  bifurcation  point  preceding  the  two  flame  fronts  was  identified  as  the 
stabilization  point.  A  Damkohler  number  based  on  the  time  scale  of  the  explosive  mode  and  the 
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scalar  dissipation  rate  was  defined  to  identify  auto-igniting  mixtures,  and  the  role  of  auto-ignition 
as  the  mechanism  for  stabilization  of  the  lifted  flame  was  discussed. 

PART  II:  DYNAMICS  OF  LAMINAR  AND  TURBULENT  FLAME  PROPAGATION 


1.  Flame  Surface  Instabilities 

We  have  extended  our  previous  work  [Jomaas  et  al.  (2007)]  on  the  diffusional-thermal  and 
hydrodynamic  instabilities  on  spherically  expanding  stretched  flames  to  include  the  study  of 
spiral  waves  which  originate  from  pulsating  instabilities  relevant  for  Le>l  mixtures  such  as  the 
rich  hydrogen-air  and  lean,  C  >  2,  hydrocarbon-air  mixtures  [Jomaas  and  Law  (2010)].  The 
experiments  were  conducted  in  our  twin-chambered  vessel  that  allows  the  spherical  flame 
propagation  at  constant  pressure.  The  experiment  involved  high-speed  imaging  at  2-25,000 
frames/s.  Five  distinct  unsteady  phenomena  were  observed,  namely,  bulk  oscillation,  target 
patterns,  ordered  and  disordered  spiral  waves,  and  extinction.  Some  of  the  key  phenomena  are 
shown  in  Fig  II- 1.  We  furthermore  have  experimentally  mapped  out  the  six  flame  responses  in 
the  phase  space  of  pressure  and  equivalence  ratio.  In  the  order  of  appearance  as  the  instability 
domain  is  approached,  entered  and  traversed,  we  observed  smooth  flames,  globally  oscillating 
flames,  flames  with  clear  target  patterns,  flames  with  mixed  patterns  that  are  intermediate 
between  targets  and  spirals,  flames  with  clear  spirals,  and  situations  of  no  ignition.  Fig.  II-2 
shows  this  map  for  H2-O2  flames.  Compared  to  FL-Air  flames  we  see  a  shift  towards  richer 
hydrogen  concentrations.  Given  that  the  difference  between  Lewis  numbers,  Le  is  not  expected 
to  be  significant  for  the  two  mixtures,  one  can  infer  that  shift  in  the  onset  of  the  instabilities  is 
caused  by  changes  in  the  reactivity  as  represented  by  the  Zel’dovich  number,  which  is  both 
directly  and  indirectly  affected  by  the  flame  temperature.  Another  interesting  observation  of  the 
H2-O2  flames  is  that  flames  at  40  and  60  atmospheres  are  actually  stable  at  very  rich  conditions. 
In  fact,  stable  flames  are  now  found  between  flames  with  instabilities  and  flames  that  extinguish. 
This  therefore  demonstrates  that  the  stability  response  can  vary  non-monotonically  with 
increasing  equivalence  ratio,  at  high  pressures.  It  is  noteworthy  that  this  observation  was 
predicted  and  explained  by  Christiansen  et  al.  (2001)  who  showed  that  the  Sivashinsky  stability 
parameter  Ze(Le- 1)  actually  first  increases  and  then  decreases  with  increasing  pressure  for 
hydrogen-air  flames.  It  order  to  evaluate  the  closeness  of  theoretical  limit  based  on  the 
Sivashinsky  criterion  with  the  experimental  results  we  have  evaluated  the  parameter  Ze(Le- 1)  as 
a  function  of  the  equivalence  ratio  for  various  pressures,  where  the  values  of  Ze  and  Le  were 
locally  extracted  from  the  response  of  unstretched  and  stretched  flames  respectively.  Fig.II-2 
then  superimposes  the  theoretical  boundary  on  the  plot  of  the  experimental  H2-O2  results,  and  it 
is  seen  that  the  experimental  transition  regime  does  span  around  the  neighborhood  of  the 
theoretical  boundary,  suggesting:  (1)  the  spiral  flame  is  a  phenomenon  involving  flame 
pulsation;  (2)  it  is  essential  to  use  realistic,  locally  extracted  thermokinetic  parameters  in 
evaluating  the  transition  criterion. 
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Figure  II- 1:  (a)  Target  patterns  on  a  spherically  expanding  lean  butane -oxygen-helium 
flame  (</>=0.59)  at  30  atmospheres.  The  physical  dimensions  of  the  frames  are  5.46cm  by 
5.46cm  and  the  photographs  were  taken  0.625  ms  apart 

(b)  Target  patterns  on  a  spherically  expanding  lean  butane-oxygen-helium  flame  ((f>=0.59) 
at  30  atmospheres.  The  frames  are  2.73cm  by  5.46cm  and  the  sequence  shows  three 
consecutive  frames  at  a  framing  rate  of  15,000 fps.  The  global  flame  appears  not  to  move  in 
comparison  with  the  rapid  changes  taking  place  on  the  flame  front. 

(c)  A  flame  photograph  of  a  hydrogen-air  flame  with  an  equivalence  ratio  of  4.30  at  20 
atmospheres,  showing  the  presence  of  spiral  waves  within  the  hydrodynamic  cells.  The 
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Figure  11-2  :  Flame  front  mode  diagram  as  a  function  of  equivalence  ratio  and  pressure  for 
hydrogen-oxygen  flames  with  the  calculated  theoretical  boundary  superimposed 


An  shown  through  many  experiments,  expanding  laminar  premixed  flames  is  generally  subjected 
to  the  excitation  of  flame  front  due  to  diffusional-thermal,  hydrodynamic,  and  buoyancy  induced 
instabilities,  which  are  respectively  caused  by  the  non-equidiffusive  nature  of  the  diffusive 
scalars  of  species  concentration  and  heat,  the  density  jump  across  the  flame,  and  body  force  in 
the  presence  of  density  gradient.  The  presence  of  cells  over  the  flame  surface  increases  its 
surface  area  and  consequently  also  the  global  propagation  speed  of  the  flame.  Furthermore,  since 
new  cells  continuously  evolve  as  the  flame  propagates,  it  is  reasonable  to  expect  that  the  flame 
speed  will  also  continuously  increase,  leading  to  the  possibility  of  self-acceleration.  In  order  to 
quantify  the  flame  acceleration  due  to  constant  evolution  of  cells  and  to  assess  the  possibility  of 
self-turbulization,  we  write  the  flame  radius  as  a  power  function  of  time,  Rav  ~  ta  ,  where  a  is  a 
constant.  Through  experiments  conducting  in  our  dual  chamber  pressure  vessel  at  high  pressure 
such  as  60  atm.  with  hydrogen  flames,  we  achieve  the  required  substantially  smaller  flame 
thickness  that  is  suitable  for  studying  cellular  propagation  [Wu  et  al.  (2013)].  Figs.  11-3 a  and  b 
show  that  experiments  with  different  ignition  energies  overlaps  with  each  other  and  the  slope  of 
the  log (Uav)  vs.  log(/?av)  remains  a  constant  until  chamber  confinement  starts  to  have  an 
influence,  at  about  Rav  =  2.0  cm.  This  result  confirms  that  a  self-similar  regime  does  exist  for 
wrinkled  laminar  flame  propagation,  after  an  initial  transition  period. 
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Rav  (cm) 


Figure  II-3:  Propagation  velocity  for  two  experiments  for  H2/air,  at  (a)  5  atm,  (/)=  0.60,  and 
(b)  10  atm,  <f>  -  LOO.  Experiments  1  and  2  are  runs  with  the  same  mixture  but  different 
ignition  energies. 
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Figure  11-4:  Experimental  acceleration  exponents  for  H 2/air  at  various  pressures  and 
equivalence  ratios  (all  data  points  are  estimated  to  have  uncertainty  +0.05;  however,  error 
bar  is  indicated  only  at  one  point  for  display  clarity) 

In  Fig.  II- 4  we  plot  the  experimentally  measured  acceleration  exponent  a  for  hydrogen/air 
mixtures  as  a  function  of  the  fuel  equivalence  ratio,  cp,  for  various  pressures.  Three  trends  are 
observed.  First,  the  values  of  a  for  lean  mixture  are  higher  than  rich  mixture  due  to  the  effect  of 
diffusional-thermal  instability,  which  is  cellularly  de-stabilizing  for  lean  mixtures  and  stabilizing 
for  rich  mixtures  leading  to  stronger  and  weaker  wrinkling  for  these  mixtures  respectively.  For 
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the  strongly  burning  mixtures,  with  (p  ranging  from  0.60  to  1.20,  we  observe  increasing  pressure 
results  in  larger  values  of  a,  due  to  reduction  in  the  flame  thickness  and  thereby  increase  in  the 
propensity  of  the  flame  to  become  hydrodynamically  unstable.  We  also  see  that  in  the  range  of  cp 
=  1.40  -  2.00,  increasing  the  pressure  does  not  change  the  value  of  a. 

We  also  find  that  for  hydrogen/air  mixtures  this  is  close  to  the  stoichiometric  value  of  <p=  1.00 
(equidiffusive  mixture),  the  acceleration  exponent  assumes  a  maximum  value  of  about  a  =  1.37, 
is  close  to,  but  still  sufficiently  smaller  than  a  =  3/2,  the  value  purportedly  related  to  the  concept 
of  self-turbulization. 

2.  Propagation  of  Turbulent  Flames 

The  turbulent  flame  speed  is  a  topic  of  wide  interest  in  combustion  and  turbulence  research  due 
to  its  importance  as  a  measure  of  flame  surface  density,  in  that  it  can  be  correlated  to  the 
volumetric  heat  release  rate  in  a  turbulent  reacting  flow.  The  problem  is  of  considerable 
fundamental  complexity,  which  is  further  compounded  by  the  disagreement  between  theories  as 
well  as  the  high  degree  of  scatter  of  the  experimental  turbulent  flame  speeds  and  their  sensitivity 
on  the  geometry  and  type  of  the  burner  used  in  the  investigation  -  a  major  hindrance  that  has 
prevented  its  utilization  as  a  meaningful  physical  quantity  for  predictions  and  for  validating 
simulations  of  turbulent  reacting  flows.  In  this  research  we  seek  a  unified  scaling  description,  at 
least  under  some  special  flow  conditions  such  as  those  in  homogenous  isotropic  turbulence.  The 
most  obvious  choice  of  the  flame  parameters  for  such  a  scaling  would  be  the  planar  laminar 

(0  /  f) 

Dcob  )  and  the  corresponding  laminar  flame  thickness^  ~  \D I cob  )  , 

assuming  a  local  laminar  flame  structure  exists,  where  D  and  are  the  characteristic  thermal 

diffusivity  and  reaction  rate  respectively  [Law  2006],  The  problem  of  turbulent  flame 
propagation  can  then  be  considered  as  a  geometric  one  in  which  the  effect  of  turbulence  is  to 
wrinkle  the  flame  at  a  multitude  of  length  scales  without  perturbing  the  inner  flame  structure. 
Such  a  problem  was  considered  analytically  in  [Chaudhuri  el  al.  (2011)],  which  showed  that  for 
a  statistically  planar  flame  propagating  in  homogenous  isotropic  turbulence,  the  turbulent  flame 
speed  normalized  by  the  corresponding  laminar  flame  speed  for  large  turbulent  Reynolds  number 
(Rer  )  is  given  to  the  leading  order  by 


sT/sL~U TTv^v^ 


ou 

l  +  \k2Y(k)dk 


1/2 


[(«„ 


(HD 


where  g  is  the  fluctuating  flame  surface  height/distance  from  the  mean  surface,  k  the 
wavenumber,  T(k)  the  flame  surface  spectrum  obtained  from  spectral  closure  of  the  G-equation 

[Peters  (1992)],  urms  the  root  mean  square  of  velocity  fluctuations  and  2,  the  velocity  integral 
length  scale  which  was  assumed  to  be  the  flame  hydrodynamic  length  scale. 
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Furthermore,  in  this  research  program  we  experimentally  measured  turbulent  flame  speeds  in 
constant-pressure  expanding  flames,  propagating  in  nearly  homogenous  isotropic  turbulence. 
Utilizing  the  self-similar  property  of  turbulent  flame  speeds  evolving  from  the  theory  and 
experimental  data  presented,  we  showed  that  the  turbulent  flame  speeds  measured  from  the 
present  spherically  expanding  flames,  as  well  as  those  from  literature  data  on  Bunsen  flames,  can 
be  scaled  by  a  single  parameter:  a  turbulent  Reynolds  number  based  on  the  geometric  and 
transport  properties  of  the  flame.  A  configurationally  independent  description  of  turbulent  flame 
propagation  in  near  homogenous  isotropic  turbulence  is  thus  proposed.  Experimentally,  we 
studied  the  flame  propagation  under  different  turbulence  intensities  for  mixtures  with  Lewis 
Number  close  to  unity  at  various  pressures.  Unlike  laminar  flames,  turbulent  flames  are  not 

symmetric  in  nature,  thus  the  meaningful  scale  is  given  by  mean  radius  (/?)  =  4  At  it  ,  where  A  is 

the  area  under  the  flame  edge  detected  by  Schlieren  images.  We  observe  the  flame  to  accelerate 
as  it  propagates  and  gets  larger,  showing  that  turbulent  flame  speed  does  not  hold  a  linear 
relationship  with  the  radius  of  the  flame. 

To  correlate  the  measured  turbulent  flame  speed  with  the  imposed  turbulence  parameters,  Eqn. 
(11-1)  provides  insight  for  self-similar  propagation  if  X,  is  replaced  with  the  relevant  largest 

flame  length  scale,  i.e.  hydrodynamic  length  scale  of  flame  surface  fluctuations  which  should  be 
a  linear  function  of  2 ft(R)  or  simply  (i?)  itself.  Consequently,  in  Fig.  II-5  we  plot  d(R)/dt, 

obtained  from  the  experimentally  obtained  (/?) ,  normalized  by  the  laminar  flame  speed  with 
respect  to  the  burnt  gas  ( SbL ),  as  a  function  of  [(unm  /  SL )((/?)/ <5^)],  i.e.  the  RHS  of  Eqn.  (II- 1 ) 

with  the  integral  length  scale  being  replaced  by  (/?) .  The  superscript  b,  as  in  SbL  implies  flame 

properties  with  respect  to  the  burnt  gas,  whereas  absence  of  a  superscript  implies  flame 
properties  with  respect  to  unburnt  gas.  Since  the  thermal  diffusivity  D  ~  SLSL ,  the  term 

\_(unm  /  SL ){{R)  /^l)]  =  Rerpe)  represents  a  turbulent  Reynolds  number  with  (/?)  being  the 

length  scale  and  the  thermal  diffusivity  replacing  the  kinematic  viscosity.  It  is  then  observed 
from  Fig.  2  that  all  the  data  from  different  conditions  of  turbulence  intensity  and  pressure,  and  at 
each  instant  of  the  propagation  event,  collapse  reasonably  well  on  a  Re“  ^  curve,  with  the 

exponent  a  =  0.54  obtained  by  nonlinear  least-square  fitting  over  the  entire  relevant  data  set. 
This  result  therefore  substantiates  the  possible  validity  of  the  1 /2-power  scaling,  as  suggested  by 
the  theory  of  [24].  More  importantly,  it  suggests  that  expanding  turbulent  flame  propagation  is 
self  similar,  at  least  in  the  domain  of  interrogation,  as  evident  from  the  relationship 

(si  r  d(R)  /  *  =  0(l)[(um,  /  SL  X(  *)  /  SL  )]'2 
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The  1/2-power  of  Eqn.  (II-2)  also  indicates  that  the  turbulent  flame  is  accelerating.  This  is  due  to 
the  fact  that  as  the  flame  expands  during  propagation,  its  smallest  wavenumber  decreases 
resulting  in  a  continuous  increase  in  the  integral  of  the  k2T  spectrum  or  the  flame  surface  scalar 
dissipation  spectrum. 


Figure  II-5:  Plot  of  flame  propagation  rate  normalized  by  laminar  burnt  flame  speed  with 
respect  to  turbulent  Re  with  average  radius  as  length  scale  and  thermal  diffusivity  as 
transport  property.  Not  included  in  the  fitting  is  the  p  /  p0=  1,  urms  ~5.57m/s  case  as  the 
Schlieren  images  show  widespread  local  extinction  in  such  flames  and  evident  from  their 
lower  propagation  rates.  ums  appearing  in  the  abscissa  is  the  u^  ((/?)),  obtained  by 

integrating  over  all  9  .  ums  -  denotes  the  average  ums  experienced  by  the  flames  during 

the  propagation  event.  The  error  bars  indicate  the  error  that  could  be  caused  by  the  mean 
flow. 


Despite  of  the  success  of 


~S~L 


u,„ 


R 


SL  SL 


1/2 


on  unity  Le  flames,  it  is  widely  recognized  that  Le 


of  most  mixtures  can  deviate  substantially  from  unity  to  affect  the  flame  temperature  and 
subsequently  the  flame  speed  through  the  intrinsic  Arrhenius  kinetics.  Therefore  it  is  of  interest 
to  seek  a  scaling  relation  which  takes  into  account  the  Le  effects  and  is  valid  for  all  fuel/oxidizer 
mixture  properties.  The  effect  of  non-unity  Le  in  modifying  the  local  flame  speed  is  most 
prominent  in  the  presence  of  local  stretch  and  curvature  due  to  the  nonequidiffusion  of  heat  and 
species  which  are  directed  normal  to  the  local  flame  surface.  The  dependence  of  the  local  flame 
speed  on  strain  and  curvature  is  essentially  nonlinear;  however,  it  has  been  shown  by  Chen 
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(2011)  and  Kelley  et  al.  (2012)  for  the  outwardly  expanding  spherical  flame,  the  empirical  linear 
relation  between  local  flame  speed  SL  and  curvature  proposed  by  Markstein  (1951)  provides  a 

reasonably  accurate  approximation,  SL  =  S,  (l  -  5 M  sr)  ,  where  k  is  the  curvature  with  the  convex 

section  in  the  direction  of  flame  propagation  being  positive,  and  8U  is  a  coefficient  having  the 
unit  of  length  and  termed  the  Markstein  length.  The  value  of  Markstein  length  in  the  unit  of 
flame  thickness  is  the  Markstein  number,  Mk  =  Su  / SL  ,  which  is  a  strong  function  of  Le.  For  a 

turbulent  premixed  flame  with  positive  Markstein  length,  the  concave  portion  of  the  flame,  with 
a  negative  curvature,  propagates  faster  than  the  convex  portion  of  the  flame  which  has  a  positive 
curvature,  leading  to  reduction  of  the  surface  fluctuations  and  thereby  the  total  flame  surface 
area.  This  dissipative  mechanism  is  most  amplified  at  large  wavenumbers  or  large  curvatures.  By 
assuming  that:  (i)  balance  of  Markstein  dissipation  by  kinematic  restoration  and  amplification  by 
thermal  expansion  due  to  their  similar  wavenumber  dependence  (but  with  opposite  signs)  [Peters 
et  al.  (2000)]',  (ii)  the  dominant  role  of  dissipation  by  Markstein  diffusion  at  large  wavenumbers 
and  thereby  retaining  Markstein  diffusion  as  the  sole  dissipation  mechanism;  and  (iii)  turbulent 
flame  speed  is  proportional  to  the  total  flame  surface  area,  the  following  scaling  relation  (Eqn. 
(69)  in  [Chaudhuri  et  al.  (2011)])  was  arrived: 
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(II-3) 


where  L,  is  the  integral  scale.  With  further  consideration  of  hydrodynamic  length  scale  for  flame 
to  be  the  brush  thickness  which  scales  with  the  mean  radius  of  the  flame,  we  arrive  at 
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Figure  II-6:  Log-log  plot  of  (a)  STx^)5  /  SL  versus  J(u^  7 SL ) ((/?)  /  <7L )  =  Re^  and  (b) 

iSV,r=o.5 1 Sl  versus  ^J(ueff  I SLh){{R) 1 8M  h)  =  Re°;5M  for  H2/air,  ^=4.0;  CH4/air  ^  =0.9; 

CoFL/air,  ^=1.3;  n-C4Hio/air  ^=0.8;  and  C2H60/air  ^=1.0  mixture  from  present 
experiments.  The  lines  represent  the  iso-CgH^  data  from  Error!  Reference  source  not 
found.. 
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Figure  II-6(a)  shows  STc=05/SL  vs.  Re^  =  {ueff  /  S^(iR)  /  SL^j  for  all  the  conditions  in  log- 
log  plots.  It  clearly  shows  that  the  normalized  flame  speed  has  strong  dependence  on  the  Lewis 


number.  Figure  II-6(b)  which  plots  STc=05/SL  vs.  Re 
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collapse  of  the  data  on  a  single  plot  proving  the  utility  of  scaling  as  proposed  in  Eq.  IF4. 


3.  Laminar  Flame  Speeds 

We  have  investigated  the  laminar  flame  speed  for  a  series  of  normal  and  cyclo-alkanes  to  build  a 
database  against  which  the  chemical  mechanism  schemes  can  be  validated.  Using  a  newly 
developed  apparatus  that  allows  experimentation  at  elevated  pressures  and  temperatures,  the 
laminar  flame  speeds  for  n-pentane,  n-hexane,  n-heptane,  and  //-octane  were  measured  at  both 
atmospheric  and  elevated  pressures  [Kelly  et  al.  (2011a)].  The  resulting  data  [Fig  77-7]  show  that 
these  n-alkanes  have  similar  flame  speeds  for  a  given  equivalence  ratio  and  pressure,  and  thereby 


18 


confirm  the  findings  of  Davis  and  Law  (1996,  1998,  1999)  while  extending  the  result  to  elevated 
pressures.  An  examination  of  the  computed  flame  structures  for  these  fuels  also  shows  similarity, 
hence  substantiating  the  suggestion  of  You  et  al.  (2009)  on  the  nature  of  such  a  similarity.  We 
have  also  demonstrated  and  explained  that  the  influence  of  stretch  decreases  with  increasing 
pressure  because  of  the  reduced  flame  thickness.  This  is  an  important  result  because  it  not  only 
indicates  that  the  systematic  error  associated  with  stretch  effects  in  the  determination  of  laminar 
flame  speeds  is  reduced  with  increasing  pressure,  it  could  also  simplify  the  modeling  and 
understanding  of  turbulent  flames  in  light  of  the  diminished  influence  of  stretch. 

Furthermore,  experimental  data  on  the  laminar  flame  speeds  of  iso-octane/air  mixtures  at 
atmospheric  and  elevated  pressures  were  acquired  using  the  counterflow  flame  and  the  outwardly 
expanding  flame,  while  the  non-premixed  ignition  temperatures  were  determined  for  an  iso¬ 
octane  pool  in  the  stagnation  flow  of  a  heated  air  jet  at  atmospheric  and  slightly  reduced/elevated 
pressures  [Kelly  et  al.  (2011b)].  These  experimental  measurements  [Fig.  II-8]  were  compared 
with  calculations  based  on  the  mechanisms  of  Curran  et  al.(2000)  and  Chaos  et  al.(2007 ),  with 
the  former  mechanism  systematically  and  substantially  reduced,  using  directed  relation  graph 
and  computational  singular  perturbation,  to  facilitate  the  calculation.  It  was  found  that  the  Curran 
mechanism  yielded  substantially  higher  laminar  flames  speeds  as  compared  to  the  present 
experimental  results  while  results  from  the  Chaos  mechanism  agree  well  with  the  present 
measurements.  These  trends  are  in  agreement  with  previous  results  on  the  laminar  flame  speeds 
for  n-heptane.  Both  mechanisms  yield  acceptable  comparison  with  the  observed  non-premixed 
stagnation  ignition  temperature. 
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Figure  II-7:  Comparison  of  flame  speeds  at  various 
pressures  with  air  as  the  oxidizer,  unburned  gas 
temperature  of  80° C.  (a)  5atm,  (b)  2  and  10  atm,  and 
(c)  Comparison  of  linearly  and  nonlinearly 
extrapolated  laminar  flame  speeds.  n-Pentane  /  air 
mixtures  at  various  pressures,  unburned  gas 
temperature  of  80° C.  Cross:  Linear  extrapolation; 
Circles:  Nonlinear  extrapolation; 
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Figure  II-8:  Laminar  flame  speed  of  iso-octane  /  air,  (a)  1  atm-298K,  (b)  latm- 
353K.  Dashed-Line:  Mech-PRF,  Solid-Line:  Mech-1 12,  Dotted-Line:  Mech-78,  (c) 
5  atm  and  (d)  10  atm. 
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